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ABSTRACT 



Using the VLA, we have performed the first observational test of dark 
matter in the form of cold, primordial fractal clouds, as envisioned by Pfenniger, 
Combes, & Martinet (1994) and Pfenniger & Combes (1994). Wc show that, 
after a Hubble Time, primordial fractal clouds will convert most of their HI 
to H2, but a small fraction of HI remains which is optically thick. This opens 
up a new window for detecting dark matter which may exist in this form. 
The detectability of such gas depends on its filling factor and temperature 
and therefore should be observable in absorption against a background source, 
with observations of sufficient sensitivity and resolution. The current VLA 
observations have made a first step towards this goal by taking advantage of a 
fortuitous alignment between the extension of the HI disk of the nearby galaxy, 
NGC 3079, and a background quasar, Q 0957+561. Our observations probe 
28 independent beams against the quasar and all of velocity space between the 
extension of a flat rotation curve and a Keplerian decline for the halo region of 
NGC 3079. We do not detect any absorption features and investigate, in detail, 
the implication of this result for the hypothesis that dark matter is in the form of 
fractal clouds. In particular, we calculate the probability that our observations 
would have detected such clouds as a function of the model parameters. The 
chance of detection is significant for an interesting region of parameter space 
(fractal dimension 1.7 < D < 2 and cloud radius 30 pc < Rc < 3kpc) and rises 
above 95% for a small region of parameter space. While our analysis does not 
rule out fractal clouds as dark matter, it does lay out the groundwork for future, 
more sensitive observations and we consider what form these might take to 
probe the range of possible cloud properties more deeply. It is interesting that 
the observations can rule out cold optically thin HI gas, if it exists, to a limit 
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of 0.001% of the dark matter. In contrast, the existence of cold HI in a fractal 
hierarchy would be an efficient way of hiding dark matter. 

Subject headings: cosmology: dark matter — galaxies: halos, individual (NGC 
3079) — quasars: absorption lines — ISM: clouds 
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1. Introduction 

The question of extended disks around spiral galaxies, heavy halos, and indeed, the 
very nature of dark matter, has engaged astronomers for many years and dates to at least 
the early observations of fiat rotation curves (e.g. Rubin, Ford, fc Thonnard 1980). Since 



that time, observations too numerous to list have been carried out in attempts to search for 
evidence, direct or indirect, for this dark matter and to constrain its nature. Whatever the 
composition, there is now strong evidence that the dark halos surrounding spiral galaxies 
can extend to large distances. According to Ashman (1992), there is general agreement 
that the dark halo of the Milky Way, for example, extends to a radius between 100 kpc and 
200 kpc. 

It has been known for some time that primordial nucleosynthesis favors a density 
in baryons that is significantly greater than the density in luminous matter (e.g., Copi, 
Schramm, & Turner 1995). While recently, some of the results of primordial nucleosynthesis 
have been called into question, the conclusion that there are nonluminous baryons in the 
Universe remains firm and is in fact supported by other arguments from astrophysics 
and cosmology (Steigman, Hata, & Felten 1999). The question remains whether baryonic 
dark matter contributes significantly to halos in spiral galaxies. Baryonic dark matter has 
attracted considerable attention of late (see Carr |1994|) , especially in the light of recent 
reports of gravitational microlensing events toward the Large and Small Magellanic clouds, 
presumably (though not conclusively) due to massive compact halo objects (MACHOs) 
in the Milky Way (Alcock et al. 1997a, 1997b). At present, microlensing experiments are 
unable to fix the MACHO fraction in the Galaxy. In addition, the most probable mass for 
the lenses (~ 0.5 — IMq) presents a puzzle since this rules out brown dwarfs and giant 
planets, initially considered to be the most plausible candidates for MACHOs. 

A number of authors have considered cold g clS clS db dark matter candidate. Pfenniger, 
Combes, & Martinet (1994) and Pfenniger & Combes (1994) (hereafter PCM and PC) have 
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argued that the flat rotation curves of spiral galaxies might be accounted for entirely by 
cold, dense H2 gas clouds. Henriksen & Widrow (1995) and Draine (1998) have suggested 
that very dense clouds might themselves act as gravitational lenses while Gerhard & Silk 
(1993, 1996) and de Paolis et al. (1995a, 1995b) have considered the possibility that gas 
might coexist with a population of MACHOs. 

PCM and PC, in particular, have proposed that dark matter consists of cold primordial 
gas clouds in extended disks or halos around spiral galaxies where the gas within an 
individual cloud is distributed in a fractal structure. Far from sources of heating, the gas 
is thermalized to the 3 K background, assisted by small quantities of H2 ice. PCM suggest 
that the slow accretion of gas toward the visible disk solves the "gas consumption problem" 
(i.e. that the timescale for gas consumption, as implied by current star formation rates, is 
shorter than the age of the galaxy) as well as the "disk-halo conspiracy" (i.e. the smooth, 
flat rotation curves through the transition region from disk to halo dominated rotation). 
Upon re-evaluation of the stellar recycling of gas, the former problem may not be as serious 
as previously thought (Kennicutt et al. 1994), but the latter problem remains. The PCM 
and PC proposal is motivated by fractal interpretations of the hierarchical structures 
observed over a wide range of scales for atomic and molecular interstellar clouds (see, 
e.g. Diamond et al. 1989, Vogelaar & Wakker 1994, and Elmegreen 1996). The HI in 7 
members of the M 81 group (Westpfahl et al. 1999) and possibly in the Small Magellanic 
Cloud (SMC) (Stanimirovic et al. 1999) have also recently been shown to have a fractal 
structure. PCM suggest that the gas is in the form of 3 K molecular hydrogen which is 
unobservable directly and, if primordial, would also be unobservable via tracers such as CO. 



However, if a component of atomic hydrogen is present (see § |2.2| ), it could be detected in 
absorption against a background continuum source with observations of sufficient resolution 
and sensitivity. 

In this paper we describe the first observational test of the PCM and PC model. 
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Specifically, we search for HI absorption associated with a foreground galaxy, NGC 3079, 
against the background quasar, Q 0957+561, which is fortuitously located nearby. Some 
preliminary results from this work have been presented in Irwin et al. (1999). In §^ we 
describe the PCM and PC model and discuss the neutral hydrogen content of the gas. In 
§1^ we describe NGC 3079/Q 0057+561, the observations and data reduction, and in §^ 
we present the results. No significant absorption features are detected. The observations 
are extremely sensitive (A^(III) = 1.4 x 10^"^ cm^) and indicate that no more than 0.001% 
of the dark matter halo can be in the form of cold diffuse HI. However, as we see in 
we are unable to rule out fractal clouds as a dark matter candidate since the probability 
of detection is high (> 95%) over a limited region of parameter space. This particular 
experiment lays the groundwork for future, more sensitive observations, as discussed in 



2. The Fractal Cloud Model 

2.1. Elementary Cloudlets 

The elementary building blocks in the model of PCM and PC are Jupiter-mass 
objects called "cloudlets" ( "clumpuscules" by PC) which describe a natural minimum 
mass for fragmentation. This scale emerges by setting the free-fall time equal to the 
Kelvin-Helmholtz time and assuming virialization. 

The cloudlet parameters, i.e. mass, radius, volume density, column density, and 
thermal line width, respectively, are given by: 

~ 4.0 X 10-3 T^/^ /i"^/^ /-^/' Mq (1) 
i?, ~ 1.5 X 10" /i-^/^ (2) 

~ 1.1 X 10^ /i^/^ / H cm-3 (3) 
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iV,~ 3.2 X 10^3 r^/4//V^' Hcm-2 (4) 
t;, ~ 9.1 X 10^2 T^/V"^^^ km s"^ (5) 

where T is the temperature, /i is the mean molecular weight, and / is a factor which 
accounts for departures from spherical symmetry and blackbody radiation (PC). Note 
that is the mass density divided by the hydrogen mass, rriH {not by ^rriH), i.e. 
ra^K = SM^/AttRIitih, and therefore includes all forms of hydrogen as well as helium. The 
same holds for the surface density, A^,.. These equations imply that warmer cloudlets are 
smaller, more massive, and thus more stellar-like. 

For T = 3K, 0.1 < / < 1 (see PC) and 1.3 < /U < 2.3 (primordial neutral gas), 
the resulting cloudlet parameters are, = 0.8-9.2 x 10"^ M©, i?* = 23-150 AU, 
= 0.25-6 X 10^ cm'^, = 0.73-2.7 x 10^^ cm'^, and = 0.10-0.14 km s-\ 



Taking n = 2.3 (neutral H2 + He, see § 2.2) and adopting / = 1, we have 



M, = 0.81 X 10"^ Mq (6) 

i?, = 23 AU (7) 

= 6.0 X 10^ H cm-3 (8) 

= 2.7 X lO^'' H cm-2 (9) 

= 0.10 km s~^ (10) 



2.2. Composition of the Cloudlets 

Cloudlets at such high densities, if they exist, should consist only of neutral gas. 
While low density galactic HI disks are expected to be truncated by the extragalactic 
ionizing radiation field (Corbelli & Salpeter 1993, Maloney 1994, Dove & Shull 1994, see 
also Bland-Hawthorn et al. 1997), a simple application of the photoionization equilibrium 
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equation shows that the background UV photon rate at the current epoch (4 x 10"^'' s~^; 
Haardt & Madau 1996) is insufficient to produce any significant ionization in such high 
density cloudlets. 

True primordial gas consists of atomic hydrogen since the "freeze-out" fraction of 
molecular hydrogen is only [H2/H] = 1.1 x 10^^ (Galli & Palla 1998) and dust is not 
available to assist in the conversion from HI to H2. However, for these high density 
cloudlets, the 3-body reactions 

H + H + H^H2 + H (11) 



H + H + H2^H2 + H2 (12) 

may be important in converting HI to H2 (see Palla et al. 1983). If we consider only these 
reactions, then the evolution of the number densities ^hi and n-a^ are governed by the 
equations: 



dt 



-27^lr^^I - 2n2ni,,nu, (13) 



dt 



Hiniii + 7^2n^InH2 (14) 



where TZi and 7^2 are the rate coefficients for reactions (0) and (0), respectively. These 
have not been measured at the temperatures and densities being considered here. Cohen 
and Westberg (1983) review the available data and recommend values, TZi = 8.8 x 10"'^^ 
cm^ s~^, independent of temperature, T, and TZ2 = 2.8 x 10^'^^ {T / K)~^'^ cm^ s~^ above 
80 K. In the absence of other data, we assume the above temperature dependencies hold 
down to 3 K so that for the case at hand, TZi is as given above and 7^2 = 1-4 x 10~^^ cm^ s^^. 
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When the cloudlets first form (time to); the gas is almost entirely in atomic form. 
However, on a timescale of order (?^hi(^o)'^i)^^ ~ lO^yrs, most of the hydrogen is converted 
to H2. Thereafter, the second reaction dominates and we have rin^ ~ ?T,Hi(io)/2 and 



For nHi(to) = 0.25 - 6 x 10^ cm"^ (§ g^) and {t - to) = 10 Gyr, Eqn. (|l|) yields nm{t) 
= 4.0 X 10^ - 9.5 X 10^ cm~^. We take nHi(t) to represent the present density of HI in 
a cloudlet, i.e. n{HI)^. This implies that the fraction of the cloudlet in the form of HI, 
f{HI) = n{HI)Jn^, is small, i.e. f{HI) = 7 x 10"^ - 2 x 10"^ 

If HI coexists with H2 throughout the cloudlet, then the HI column density will be 
N{HI)^ = 2r2Hi(t)i?* = 2.8 x 10^^ - 4.3 x lO^o cm'^. Since in general. 



where r„ is the HI optical depth at velocity, v, and T, is the spin temperature (kinetic 
temperature) of the gas, the above parameters imply that r ~ 5 — 700, and hence even 
this small fraction of HI in the cloudlets would be optically thick. The detectability of such 
gas depends not on quantity, but on temperature and filling factor. Thus, even though HI 
may constitute a small fraction of the mass of the cloudlet, it can still provide a sensitive 
observational probe of the gas, if it exists. 



nniit) 



1 



(15) 



'J^2nm{to) (t-to) 




(16) 



2.3. 



Fractal Clouds 



Cloudlets coalesce hierarchically in a fractal distribution to form larger structures 
called clouds (denoted with subscript, c, below). The clouds have parameters related to the 
parameters of the cloudlets via. 
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M, = KM, = (_£j M, (17) 



D-3 



Re 

nr = \ — \ 



= (^-^j iV, (19) 

= (-^j V. (20) 

where A/'c is the number of cloudlets in a cloud, Vc is the velocity dispersion of the cloud, 
Rc/ -R* is the dynamic range of the fractal structure, and D is the fractal dimension. A pure 
diffuse smooth medium has D = 3 and lower values of D correspond to progressively lower 
filling factors. In general, we expect D ^ 2 since for D > 2, the clouds tend to dissolve 
due to internal collisions (PC). Molecular clouds in the ISM exhibit fractal structure with 
1.4 < D < 2. Recent results for HI in the M 81 group of galaxies give 1.2 < D < 1.5 
(Westpfahl et al. 1999) and for HI in the SMC, D = l.h (Stanimirovic et al. 1999). 



3. Observations and Data Reduction 

3.1. The NGC 3079 - Q 0957+561 Pair 

A fortuitous alignment of the extension of the major axis of the galaxy, NGC 3079, 
with the background quasar, Q0957+561 has provided the opportunity to search for cold 
gas in an extended halo or an extended disk around the foreground galaxy (see Figure 1). 
NGC 3079 is at a distance of D = 15.6 Mpc (Hq = 75 km s^^ Mpc"^; z = 0.004) and has 
been mapped in HI by Irwin & Seaquist (1991) and Irwin et al. (1987). The outermost 
point at which emission is observed is denoted by a small star in Fig. 1. Q0957+561 is 
a gravitationally lensed system ( Walsh et al. 1979|) at a redshift of z = 1.4 and has also 



been extensively studied (e.g. Schmidt & Wambsganss 1998 and references therein). The 
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projected separation between the center of NGC 3079 (taken to be the kinematic center 
given by Irwin & Seaquist 1991), and Q0957+561 (taken to be the position of component 
A of RA (J2000) = 10^ Ol"' 20?69, DEC (J2000) = 55° 53' 55!'9) is 14.13 arcminutes, or 
64.1 kpc at the distance of NGC 3079. For our purposes, we consider the lensed quasar to 
be simply a source of background radio continuum emission. 

3.2. The Observations 

Observations were carried out in the 21 cm spectral line on 03 January 1997 in the 
A configuration of the Very Large Array (VLA)|[]. The flux calibrators were J1331+305 
(3C286) and J0137+331 (3C48) which had flux densities at the central observing frequency 
of 14.79 Jy and 16.04 Jy, respectively, using the latest VLA calibration scale. The phase 
calibrator, J0957+553, for which we find a flux density of 2.45 ± 0.01 Jy, is separated from 
the source by only 43 arcminutes and was observed every ~ 20 minutes. All calibrators 
had UV restrictions which were applied during calibration. The field center was placed 13 
arcseconds south of the quasar to avoid baseline-dependent map errors which tend to show 
up at the field center (Ekers 1989) and, for highest sensitivity, the 4IF mode was used. 

The systemic velocity of NGC 3079 is 1117 km s~^ (heliocentric, optical definition) and 
the peak of the rotation curve is offset ±215 from this value (Irwin & Seaquist 1991) with 
the background quasar occurring on the blueshifted side of the galaxy. Thus the central 
observing velocity was set to 974 km s~^, which is between the velocity of a flat rotation 
curve for NGC 3079 extended to the projected distance of the quasar (i.e. 902 km s~^) 
and the velocity expected at the position of the quasar for a Keplerian fall-off (1010 km 
s^^). The bandwidth of 163.4 km s~^ encompasses both of these possibilities (though some 

^The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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channels at the band edges were lost due to a loss of sensitivity). At these frequencies, there 
should be no contamination by Galactic HI or by the lensing galaxy (z = 0.36) or known 
absorption systems (z = 1.39 and z = 1.12, Walsh et al. 1979| ) along the line of sight. 



With on-line Hanning smoothing and 63 channels, the final channel width (= the 
spectral resolution) was 2.593 km s~^. The primary flux density calibrators were also used 
as bandpass calibrators and were observed 4 times during the 10 hours of observations. The 
calibration observation closest in time to the source observation was applied. Examination 
of the average bandpass for each time showed that the maximum temporal change in the 
bandpass over the 4 scans was of order 0.2%. 

UV data for each channel were Fourier Transformed and cleaned (e.g. Cornwell & 
Braun 1989) using uniform weighting (e.g. Sramek & Schwab 1989) (RA-DEC- Velocity) 
cube. A uniform weight clean continuum image was also made by averaging together 
the UV data for all channels (except the noisy ones at the band edge) before Fourier 
Transformation. The data were phase-only and then phase-|-amplitude self-calibrated (e.g. 
Cornwell & Fomalont 1989) until no improvement in the rms resulted. The final rms noise 
per channel is 0.61 ± 0.01 mJy beam"^, where the error represents the la variation over all 
channels. This is ~ 1.3 times the theoretical rms for this weighting. 

A continuum-subtracted cube was made by subtracting a linear fit to the visibilities 
over low-noise channels from each total emission channel. The resulting rms noise in the 
continuum-subtracted maps is 0.58 ± 0.01 mJy beam"^. Channel-to-channel variations due 
to the bandpass calibration are smaller than this. The spectral dynamic range, i.e. the 
ratio of peak continuum intensity to rms noise in the channel maps is ~ 300/1. 

Naturally weighted cubes and continuum images were also made (with and without a 
UV taper), but since the rms was not significantly improved and the beam size was larger, 
subsequent analysis proceeded on the uniform weighting maps. 
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A summary of the observing and map parameters is given in Table 1. 



4. Results 

The continuum map is shown in Figure 2, with components labelled according to 
Greenfield et al. (1985) and Avruch et al. (1997). The lensing galaxy is associated with 
radio source, G, which is visible as a northern extension to component B. This map is very 
similar to the 18 cm image obtained by Avruch et al. (1997) from archival VLA data and 
confirms their discovery of sources Rl, R2, and the northern source, N. Although not visible 
in Fig. 2, we also detect their southern source, S, in the naturally weighted map. 

The continuum-subtracted cube was inspected for signs of emission or absorption 
features. No HI features were detected. In Fig. 3, we show spectra taken at the positions of 
the peaks of continuum components A, B, C, D, and E. The last panel shows the total fiux 
as a function of velocity, integrated over a region in which the continuum emission is greater 
than 10 X the continuum map rms noise value. These show no evidence for features above 
~ 3 X the rms noise level per channel. Note that the noise in these spectra is dominated by 
the rms map noise, rather than by variations in the bandpass calibration which are smaller 
(§ 3^). There is a slight hint of an absorption feature near zero velocity for components A, 



C, and E. However, this feature disappears if a spectrum is taken several pixels away from 
the continuum peaks but still within the same beam, and it also disappears if an average 
spectrum is taken over the beam at these positions. Moreover, if a spectrum is taken at 
a random position off-source, similarly sized features result. Thus, these features are only 
noise peaks. 

We also formed 0th and 1st moment maps, i.e. total intensity and intensity-weighted 
mean velocity maps, respectively (not shown) in which low intensity emission was cut-off 
so that the higher intensity peaks would be highlighted. Again, there is no evidence for 
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emission or absorption in any channel. 

In the event that HI exists at a low level over many channels, we also created a 
continuum-subtracted cube by doing linear fits to the visibilities over channels near the 
(low-noise) ends of the useable band, only. No evidence for any HI features could be seen 
in this cube either. 



5. Constraints on Fractal HI in a Dark Matter Halo 

5.1. Density of HI within NGC 3079 

Our abihty to detect clouds of the type proposed by PCM and PC depends both on 
the structure of individual clouds (the parameters R^, i?*, and D) and the overall mass 
density in clouds. To quantify the latter, consider the hydrogen clouds in NGC 3079 to 
be smoothed into a continuous mass distribution throughout the halo of the galaxy. The 
line of sight to the background QSO intersects various galactocentric radii of this mass 
distribution. We define the hydrogen mass density at the radius closest to the center of 
NGC 3079 along this line of sight to be Psystem- 

Flat rotation curves observed in the outer parts of spiral galaxies suggest that dark 
matter halos follow an density law. Beyond this, little is known about the distribution 
of dark matter (e.g., shape of the halo). For definiteness, we assume an isothermal spheroid 
model for the halo of NGC 3079: 



~ AnG + zyq^ ^^^^ 

where Vc is the circular rotation speed, R is the galactocentric radius in the plane of the 
disk, z is the coordinate perpendicular to the plane (r^ — R^ + z"^), q is the flattening 
parameter (g < 1 for an oblate halo), and A(g) = y/1 — j {q arccos q) is a geometric factor 
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equal to 1 for q = I and rising to 6.8 for q = 0.1. We have ignored the possibihty of a core 
radius since we are only interested in the outer parts of the halo. 

Roughly speaking, a maximum value pmax for Psystem is obtained by assuming that the 
entire halo is composed of clouds (i.e., Psystem = Pdm), that the halo is highly flattened 
(g = 0.1) and that the projected position of the quasar is in the equatorial plane of the 
halo {R = 64kpc; z = 0). Taking Vc = 215 km s^^ (Irwin & Seaquist 1991) we find 
Pmax = 9.6 X lO^^^gcm"^. Let = Psystem/Pmax bc the mean density in clouds relative to 
Pmax- Fh = 1 therefore corresponds to the "best-case" scenario for detecting clouds. Any 
possible effect which might reduce Psystem is parametrized by Fh- 

Integrating Eq. ( pi]) from -oo to oo along the line of sight to the QSO and expressing 
the result in terms of a maximum surface density A^max (i-e., the surface density assuming 
g = 0.1 and z = 0) and the parameter, Fh, yields the following for the surface density: 



A^system = N^.. Fh = J^^^H = 3.6 X 10^' Fh Hcm^^ (22) 

The units, Hcm~^, are the same as in Eq. (^ and therefore represent the mass per unit 
surface area normalized by mn, rather than the number of molecules per unit surface 
area. In deriving this expression, we have assumed that the density law is well beyond 
R = 64kpc. Obviously, one could introduce an additional parameter to allow for the 
possibility that the density law is steeper than in the outer parts of the halo. 



5.2. Observational Parameters and Their Relation to Fractal Parameters 

5.2.1. The Observational Parameters 



For HI absorbing gas in front of a uniform background continuum source, the measured 
brightness temperature, Tg, in a given beam and velocity channel (i.e. a given resolution 
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element) is given by 



Tb = [Tce-^ + T,(l - e-^)\ MhU + T,{1 - N h U) (23) 

where Tc is the brightness temperature of the background continuum, Tg, t, are the spin 
temperature and optical depth, respectively, of the gas, fh, are the filling factors for 
the beam area and velocity for a single cloud, and A'' is the number of clouds in the 
resolution element. Thus, the quantity, N" fbfv is the filling factor for a resolution element. 
The factors, fi, and fy, are dealt with separately since there are different conditions for 
saturating the beam and channel, respecively. The above equation reduces to the usual 
one for HI if this filling factor is unity. Aside from a constant of proportionality given by 
the Rayleigh- Jeans relation (T/I, Table 1), Tb is represented by the total emission channel 
maps. 

After continuum subtraction, 

ATb = (T, - T,)(l - e--)Nhfv (24) 

where ATb represents the measured quantity in the continuum-subtracted cubes. 
Since for most emission, Tc » Ts, then we can form the ratio. 



ATb ^ fv fb foi' T > 1 (fractal clouds) 

rp \ (^^) 

y T fv fb for r < 1 (optically thin gas) 

Figure 4 shows a map of the quantity a = ATb/{—Tc) for a single velocity channel where 
we have taken ATb to be the negative of the rms map noise. To ensure that T,. >> T^ 
for all points, we have cut off all continuum emission less than 1.5 mJy beam~^ = 583 K. 
Since the rms noise is constant over the map, this map has the appearance of the reciprocal 
of the continuum map. Such a map could be formed for any velocity channel, but the map 
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for any other channel would be identical, since both the continuum distribution and rms 
noise are the same, channel to channel. Note that this ratio map would be a map of upper 
limits to the optical depth under the more common assumptions of unity filling factor and 
optically thin gas. 

Figure 4 provides observational constraints on the filling factor of optically thick fractal 
clouds through Eqn. (]25|). If a is the value read from Figure 4 for a particular resolution 
element, then our non-detection implies 

Mhh<?>a (26) 

where we impose a 3cr upper limit, i.e. [3 a], on detectability. 

The minimum of Figure 4 has a value 0.0033. Thus the strongest constraint we have 
for a single beam is 3a = 0.01, i.e. less than 1% of this beam/velocity resolution element 
is filled with cold, fractal clouds. Of course, the clouds are distributed randomly and so 
we can tighten our constraint by considering the full data set. Our observations cover 
approximately 28 independent beams in any given channel and there are 163 channel^ or 
4564 independent beam/velocity resolution elements, in total. 

5.2.2. Model Parameters 

In this section, we apply the observational constraints to a family of fractal models 
defined by the parameters Rc/R^, D, and Fh. We assume that the cloudlets are optically 
thick in HI (§ |2.2| ) and that their physical characteristics are given by Eqs. (|6|)-(|10D. We 
focus attention on values for D between 1 and 2 ([ ^2.3| ) and values for Rc/R* less than 

^We have taken the full velocity coverage in this analysis, rather than the slightly reduced 
coverage due to the presence of a few noisy end channels. 
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3 X 10^ which corresponds to a cloud radius ~ 3kpc. It is difficult to imagine a cloud 64kpc 
from the center of the galaxy with a size much bigger than this. Indeed, such a cloud would 
subtend 2.7° if placed around the Milky Way, albeit with low covering factor. 

For each point in parameter space, we calculate the probability that a cloud would have 
been detected. Recall that for a positive detection in a given resolution element, we require 

fbfv > 3a. It is straightforward to determine fb and as a function of Rc/R*, D, 
and Fh (see below) while N can only be determined in a statistical sense. To see how 
this works, consider first a single beam. Let i-*(Abeam) be the probability of finding Abeam 
clouds in this beam and -P(A/'beam|.A/'eff) be the conditional probabihty that at least one of 
its channels contain Nes or more clouds, given A/beam- By Bayes's theorem, the probability 
that at least one channel will contain HeS or more clouds is 



P{Ks) = P(A'beam)P(A'beam|A;ff) (27) 

The probability of having a detection in this beam is therefore 



P= E P{K^) (28) 

A/'eff=A/'min 



where 



A/min = 4^ (29) 

JbJv 

is the minimum number of clouds required for a detection. 

It is straightforward to extend this analysis to all 28 beams. We use j =1-28 to label 
the different beams, i.e., P, a, and Hrain all carry the subscript j. The probability for 



- 19 - 



detecting at least one cloud within the entire data set is 

V = l-lli^-P,) (30) 

j 

The calculation of V requires, for each point in parameter space, f^, fb, -PlAbeam), and 
-P(A/boam|A/'efr). We determine these now, except for P(A/boam|A/'efr) which is presented in the 
appendix. 

Recall that the velocity dispersion of a cloud is Vc given by Eq. (|2|) while the width of 
a velocity channel is Av = 2.59 kms~^. A contour plot of Vc as a function of Rc/R* and D 
is given in Figure 5. We estimate by assuming that for > Av, the cloud covers the 
entire channel while for Vc < Av the cloud covers a fraction Vc/Av of the channel: 

fy = min{l, vJAv} (31) 

This ignores the obvious complication that a cloud can straddle more than one velocity 
channel. In particular, for Vc — Av, a cloud will generally cover some fraction of two 
neighboring channels. 

The beam filling factor, depends on both the size of the cloud relative to the beam, 
and the area filling factor for the cloudlets in a cloud. If the cloud is larger than the beam 
{Rc > Rb where Rb = 47 pc is the mean beam radius) then fh = {Ri,/ R^)^'"^ , in analogy 
with Eq. (ID. Conversely, if R^ < Rb then fb = (AJAb) (NJN^) = {RjRbf (Rc/R*)^'^ 
where Ac and Ab are the areas of the cloud and beam respectively. These results can be 
summarized in the following expression: 



where -Rmin = min (i?b, Rc). From Eqs. (p9D,(|3TD and (|3^) we construct a contour plot 
(Figure 6) for Afmm assuming a = 0.0033 (our most sensitive beam). From this figure, 
we can read off the number of clouds required in our most sensitive beam for a positive 
detection. For example, if the cloud system is described by the parameters D = 1.64 and 
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Rc/R* = 10^, then we would require at least one cloud in this beam for a positive detection. 
Likewise, for D = 1.47 and Rc/R* = 10^, we would require about 10 clouds in this one 
beam. 

We now calculate (Abeam), the expectation value for the number of clouds in a single 
beam. To do this, consider an area A large enough to contain a representative sample 
of clouds. The expected number of clouds in A will be Ma = (Nsystem/Nc) {A/ Ac). Next 
consider a single beam in A. If R^ > Rb there are three possibilities: cloud and beam do 
not intersect; cloud partially covers beam; cloud covers the entire beam. Likewise, the three 
possibilities for Rc < Rb are: cloud and beam do not intersect; part of cloud covers beam; 
cloud entirely within beam. To avoid a rather awkward treatment that would attempt to 
take into account each of these possibilities, we use, as criteria for the cloud to be "in" the 
beam, that either the center of the beam is within the cloud perimeter (the Rc > Rb case) or 
the center of the cloud is within the beam perimeter (the Rc < Rb case) . The probability for 
a single cloud in A to be also in the beam is therefore A^^^^/A where A^ax = max(Ab, Ac). 
The probability that none of the clouds in A overlap with the beam is (1 — At^s.x/ A)"^-^ ^ 
(1 — Ma Aiaayi/A). Since A > Ma ^max, the expected number of clouds in the beam is then. 



/ \r \ \r A I A -^system /^max\ -^H-^max f Rc\ f Rma,x \ /oo\ 
(A/beam) ^ Ma A^^^/A = — — = — — — (33) 



where -Rmax = niax {Rc,R^). 

A contour plot of (Abeam) is given in Figure 7. We have assumed that Fh = 1, but the 
contours can be shifted proportionately for lower Fh according to Eqn (^). Note, however, 
that Fh approaching 1 is the only astrophysically interesting case. 

The desired probability function, P(A/beam) is given by the Poisson distribution: 

^(A/'beam) = ^ '"^^V ■ ^^4) 

-' Vbeam • 

The calculation of -P (A/beam I A/'eg) is more involved and therefore left to the appendix. 
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All of the ingredients are now in place and we can perform the calculation for V, the 
probability for detecting at least one cloud in the entire dataset. Recall that an oblate 
halo filled with clouds, i.e. Fh = 1, is the best case scenario for detection. This is shown 
in Figure 8. The consequences of the beam area filling factor, Eq. (|3^), are evident in the 
plot. For example, clouds with D < 1.7 have < 0.02 even if they fill the beam {Rc > Rb)- 
Likewise, clouds with Rc < O.I-R5 ^ lO'^'^i?* have < 0.01 even for D = 2. In these regions 
of parameter space one requires two or more clouds in a given resolution element which is 
unlikely. If fractal clouds exist in this range of parameter space, our observations would not 
have been able to detect them. 

Our strongest constraints apply to the region, 1.75 < D < 1.85 and 100 pc < Rc ^ 3kpc 
{Rc > Rb) for which we find the highest probability, i.e. ~ 95%. Thus, it is likely that if 
such fractal clouds exist, they would have been detected in our observations. A small cloud 
in this region, e.g. D = 1.8, -Rc = HO pc, would have a mass, 5 x 10"^ Mq (Eqn (|1^) 
and a velocity dispersion of 25 km s~^ (Eqn (|20|)). There is another interesting region 
of parameter space which can be ruled out at the 90% probablity level. This extends to 
much smaller clouds {Rc < Rb) with a higher covering factor (higher D). For example, 
D = 1.9, Rc = 17 pc (i.e. log(-Rc/-R*)=5.2) correponds to a 6 x 10^ Mq cloud with a 
velocity dispersion of 22 km s~^. Over all, there is a region of parameter space with a fractal 
dimension between about 1.7 and 2 and a fractal cloud dynamic range between about 10^ 
and 10^ within which there is a greater than 50% probability that a detection would have 
resulted from our observations. 

5.2.3. Comparison with Optically Thin, Diffuse Gas 

It is fruitful to examine the more conventional assumption that the gas is optically thin, 
diffuse, and that the velocity and beam filling factors are both 1. In this case. Fig. 4 provides 
upper limits to r. Our strongest constraint, i.e. r < 0.01, then gives, for a single channel 
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(Eqn dg)), A^(ifJ) = 4.7 x lO^^ T, cm-^ and for T, = 3K, N{HI) = 1.4 x 10^^ cm.-^, 
illustrating the high sensitivity of these absorption observations. This limit is deeper than 
previously measured 3a noise values of N{HI) = 1 x 10^^ [Tg/ fb)cm~^ for Lya absorption 
systems (Carilli et al. 1998), and is comparable to the upper limit of 3 x 10^^ Tg/fb for 
HI absorption against more distant cooling flow galaxies (Dwarakanath et al. 1994). From 
Eqn (|22D, this implies that Fh = 6.4 x 10^^ for a flattened halo (and a factor of 7 higher 
for a spherical halo). Thus no more than 0.001% of the dark matter halo in NGC 3079 can 
be in the form of cold diffuse HI. 

6. Discussion 

The detection of dark matter remains one of the most ambitious quests in modern 
astrophysics. The observations described above demonstrate that one can search for a 
particular type of dark matter, cold, dense, hydrogen clouds, by looking for absorption 
features in the continuum spectrum of a background quasar. This is, of course, reminiscent 
of the rich study of Lyman-a clouds. However, in the case of NGC 3079/Q 0957+561, there 
is no doubt that the foreground absorber is a galaxy, in this case one which is nearby and 
has been well-studied in emission. We are also probing, in particular, the possible existence 
of halo gas far beyond the emitting HI disk which may be in an unusual form, i.e. cold and 
with a fractal structure. 

The analysis indicates that for an all-cloud halo with model parameters in a limited 
though interesting range (fractal dimension D > 1.7; dynamic range of fractal ~ 10^ — 10'') 
there was a good chance for detection (see Fig. 8). The areas for improvement are obvious. 
Improved resolution would extend the region of sensitivity to lower values of Rd R^,- Better 
signal to noise would bring probability levels up to where one would be able to rule out, 
at say the 99% confidence level, certain regions of parameter space. In addition, improved 
signal to noise would extend the region of sensitivity to lower values of the fractal dimension 
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D. Roughly speaking, we require {Rb/ R^)^^"^ > 3a where a = (Tb — Tc)/{—Tc). An 
improvement in signal to noise by a factor of 10, at the same resolution, would allow us to 
probe clouds with D as low as 1.45. 

Since these observations were at the highest resolution possible for connected arrays, 
the only practical improvements would be to search for objects like this in our own Galaxy, 
to use VLBI techniques, and/or to argue for more observing time for higher S/N. The 
difficulty in probing for objects in our own Galaxy is that the distance to the absorber may 
not be well known so it may not be possible to pin-point a particular absorber as a halo 
object. In moving to VLBI techniques, sensitivity is lost in comparison to connected arrays. 

One possibility might be use VLBI techniques to look for individual Jupiter-mass 
cloudlets, directly, in our own Milky Way halo. Consider, for example, a small "sparse" 
cloud with D = 1.3 and {Rc/R^) = 10^ (-R* = 23 au, Rc = 0.1 pc), at a distance of 
30 kpc which is distant enough to be far from Galactic sources of heating. An individual 
cloudlet will subtend an angle of 1.5 milliarcseconds while the cloud will subtend an angle, 
1.5 arcseconds. The mean distance between clouds will be [r,^ /M^ ^ = 1 x 10^^ pc (see 
Eqn (|13)); or 8 milliarcseconds. These parameters are favourable for multiple detections 
over a VLBI field of view against a strong background source. If the cloudlets are indeed 
optically thick at 3 K, the absorption signature should be obvious. 

Indeed, VLBI HI absorption measurements have already detected a population of 
objects in the Galaxy with size scales of order 10s of au, as summarized and discussed by 
Heiles (1997). These tiny HI clouds, referred to by Heiles as "tiny scale atomic structure" 
(TSAS) have densities of ~ 10^ cm~'^. This is in the range of 10^ 10^ cm~^ which we 
have estimated for the HI fraction of PC cloudlets ( ^2.2| ). TSAS features are relatively 
nearby as indicated by velocities, \v\ < 20 km s^^ (Diamond et al. 1989, Davis et al. 1996) 
or by absorption against pulsars which are within ~ 3 kpc (Frail et al. 1994). Consequently, 
they are near Galactic sources of heating (temperatures are often assumed to be ~ 50 i^). 
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Given the sizes and densities of these objects, the possibihty is raised that they are PC 
cloudlets which have already been accreted onto the Milky Way. A clear detection of CO 
associated with the HI absorbing systems would rule out this possibihty. 

To resolve the problem of unknown distance and argue that such objects are in the 
Galactic halo, a population of such objects would have to be detected at high velocities 
over a number of lines of sight. 

7. Summary 

Pfennigcr ct al. (1994) and Pfenniger & Combes (1994) have proposed a new model 
of dark matter consisting of cold clouds, radiatively coupled to the 3 K background, in 
an extended disk around galaxies. This gas is proposed to consist of virialized elementary 
cloudlets and to be of primordial composition except that at their high densities (~ 6 x 10^ 
cm~^), gas phase reactions can convert most of the HI to H2. The cloudlets are arranged 
in a hierarchical fractal structure into larger clouds, the fractal structure being essential 
to ensuring that collisions do not dissipate the cloud and, in combination with low 
temperatures, that stars do not form. 

In this work, we have shown that, at the temperatures and densities envisioned, a small 
fraction of HI will remain in the cloudlets and will likely be optically thick. Thus, their 
detectability depends on their temperature and the filling factor of the resolution element, 
rather than quantity. At a temperature of 3 K, the clouds cannot be seen in emission 
against the microwave background, but HI is potentially observable in absorption against a 
background source. This means that HI can be used as a probe of PC-type fractal clouds 
and opens a new window for exploring halo dark matter which may be in this form. 

Using the VLA, we have carried out the first observational test of the PC model by 
taking advantage of a fortuitous alignment of a strong background quasar at a projected 
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separation of 64 kpc from the center of the nearby galaxy, NGC 3079. The quasar is 
ahgned along an extension of the galaxy's major axis. By tuning to the velocity of the 
foreground galaxy, we made 21 cm absorption measurements against the quasar, obtaining 
high spatial resolution (1.2", or 28 beams across the background continuum) and high 
velocity resolution (2.6 km s~^). We do not detect any HI absorption to a (minimum) 3cr 
upper limit 3ATbJ{-T^) = 0.01. 

Our observational limits can be related to the filling factor of fractal clouds in a 
beam/velocity resolution element and we have considered, in detail, how the observations 
might constrain parameter space, defined by the quantities, Rc/R*, D, and Fh, i.e. the 
dynamic range of the fractal structure, the fractal dimension, and the fraction of dark 
matter contained in such clouds, respectively. We find that much of parameter space 
cannot be ruled out by our observations. However, there is a region of parameter space 
(i.e., 1.7 < D < 2 and 30 pc < Rc ^ 3 kpc) within which there is a reasonable probability 
of detection. In the region 1.75 < D < 1.85 and 100 pc ^ Rc ^ 3 kpc, the probability of 
detection is as high as ~ 95% (see Fig. 8). Thus, it is unhkely that clouds with these 
parameters exist. 

It is interesting that cold diffuse (optically thin) HI can be ruled out to a hmit of 
0.001% of the dark matter. In contrast, congregating the HI into optically thick fractal 
clouds is a very efficient way of hiding dark matter. Our analysis, applicable to any fractal 
clouds containing optically thick HI, is the first step in probing dark matter in this form. 
We outline possible ways of improving the observations in the future. 

The authors wish to thank T. Abel, M. Mandy, R. Henriksen, N. Pelavas, and M. 
Walker for useful discussions. We have made used of the NASA/IPAC Extragalactic 
Database (NED). This work is supported by the Natural Sciences and Engineering Research 
Council of Canada 
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Appendix 

In this appendix we calculate P ( Abeam |A^ff), the probability that at least one channel 
in a given beam will contain jVcs or more clouds given jVbeam clouds in this beam. The 
clouds are characterized by an internal velocity dispersion, Vc (assumed to be the same 
for all clouds) and a bulk velocity, v which is different for each of the clouds. For a given 
channel to contain contributions from two clouds, we require that their bulk velocities 
be separated by less than max(t'c — Av,Av). That is, for Vc > 2Av, we require that 
the overlap of the clouds in velocity space should be greater than the width of a single 
channel whereas for Vc < 2Av we require that the separation of the clouds in velocity space 
be less than Av. To continue, we need to know the distribution function for v. In the 
absence of a detailed model for the cloud population, we assume that the v are equally 
likely to reside in any part of the 163 km velocity range covered by our observations. 
The probabihty that two particular clouds will contribute to a single channel is therefore 
p = max(vc — Av, At')/163kms~^. Put another way, we can think of there being L = 1/p 
bins in velocity space and we are asking for the probability that a certain number of clouds 
will appear in a single bin. 

In general, if there are Abeam clouds in the beam, then the the probability that Iq bins 
contain clouds, li bins contain one cloud, etc is given by 



P(A-beam, p) = P^m ^^j^;: (M) 

i-i-ml'm- {ml) 

The first factor is simply the probability of having a particular combination of Afbeam clouds 
in L bins. The second factor handles the fact that we don't care which bins contain the 
different numbers of clouds. The final factor handles the fact that we don't care which 
clouds go into which bin. The desired probability is then given by 



(A2) 



where the sum is over all combinations of Im consistent with the constraints that 

Em Im = Em ^m"^ = A/'beam and > 1- 

To illustrate how the calculation will go, we consider a few simple cases. Clearly 
-P(Abeam|l) = 1 since the probability of having at least one channel with one or more 
clouds is unity provided Abeam > 1. As discussed above, and in agreement with EqAl 
and A2, P(2|2) — p. For P(4|2), there will be four contributions to Eq. A2 corresponding 
to {lo, h, h, k, k) = (1, 2, 1, 0, 0), (2, 1, 0, 1, 0), (3, 0, 0, 0, 1), (2, 0, 2, 0, 0). For the special 
case where I = 1/p — A we find -P(4|2) = 29/32. This case is easy to check since 
the probability of having two or more clouds in one of the four channels is simply 
1 — (probability of having one cloud in each channel) = 1 — 4!/4^ = 29/32. 
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Fig. 1. — Digitized Sky Survey optical image of the foreground galaxy, NGC 3079, and 
the background quasar, Q 0957+561 (large star). The small star marks the approximate 
boundary of the HI emission associated with NGC 3079 as measured by Irwin et al. (1987). 

Fig. 2. — Uniformly weighted continuum map of Q0957+561. Contours are at -0.26, 0.26 
(2(t), 0.3, 0.4, 0.5, 0.8, 1.5, 2.5, 5, 10, 20, 30, 50, 75, 120, and 170 mJy beam-^. The peak 
intensity is 174.5 mJy beam~^ and the beam size is 1.31" x 1.15" at position angle, -87.31°. 
Labelling follows the conventions of Greenfield et al. (1985) and Avruch et al. (1997). The 
lensing galaxy is associated with radio source, G. 

Fig. 3. — Spectra through the continuum-subtracted channel maps at positions 
corresponding to the peak of components A, B, C, D, and E (see Fig. 2). The rms map 
noise per channel is 0.58 ± 0.01 mJy beam~^ and the velocity is with respect to 974 km 
s"^. The last panel shows the integrated flux (mJy) measured over a region within which 
the continuum emission is greater than lOcr. 

Fig. 4. — Map of upper limits on fv (if the gas is optically thick) or on r (if the gas 
is optically thin) for a single channel, assuming that the upper limit on ATb is given by 1 
X the rms map noise. The map values range from a minimum of 0.00332 to a maximum of 
0.385 and the greyscale ranges from 0.003 (black) to 0.06 (white). External black areas have 
been blanked. 

Fig. 5. — Contour plot of the cloud velocity dispersion Vc (in units of kms"^) as a function 
of model parameters Rc/R* and D. Contours are as labelled. The 2.59 and 163 contours 
correspond to Vc — Av and Vc equal to the velocity range of the observations respectively. 

Fig. 6. — Contour plot of A/'min- Contours are at equally spaced intervals in /ogiio(A/'min)- The 
horizontal line corresponds to the beam radius. 
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Fig. 7. — Contour plot of (A/beam)- Contours are at equally spaced intervals in /ogfio ((Abeam) )■ 
We have assumed Fh — I. 

Fig. 8. — Contour plot of the probability V for detecting a cloud as a function of the model 
parameters Rc/R* and D. For this figure, we assume Fh — 1. Probability contours are as 
labelled. 
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Velocity (km/s) 
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Table 1. Observing and Map Parameters 



Observing Date 
On-source Observing time 

^center {ycenter^ 

Av {\v) 

No. of spectral channels 
Bandwidth 

Synthesized Beam Size 

Rms Map Noise/channel [Total emission] 

Rms Map Noise/channel [Continuum-Subtracted] 

Rms Map Noise [Continuum map] 

T/I 



03 Jan 1997 

gh 54m 

974 km s-i (1415.8563 MHz) 
2.593 km s"! (12.207 kHz) 
63 

163.4 km s-i (769.0 kHz) 

1.31" X 1.15" @ -87.3° 
0.61 ± 0.01 mJy beam-i 
0.58 ± 0.01 mJy beam"! 
0.13 mJy beam~^ 
404 K/(m.Ty Ijeam-i) 



Notes to Table 1. 

All velocities heliocentric, optical definition, ^center is the center of the observing 
band. Av gives the channel separation = the velocity resolution. 



